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t>" • Abstract 

o 

CN| I Motivated by the recent measurement of low longitudinal polarization fraction in the 

decay mode B® — ► <pK*°, which appears not to be in agreement with the standard model 
expectation, we analyze this mode in the minimal supersymmetric standard model 
with mass insertion approximation. Within the standard model, with factorization 

M-t, 

approximation, the longitudinal polarization is expected to be Jl ~ 1 — 0(l/ml). 
We find that this anomaly can be explained in the minimal supersymmetric standard 
model with either LR or RL mass insertion approximation. 
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1 Introduction 



One of the important goals of the B-factories is to verify the standard model (SM) predictions 
and to serve as a potential avenue to reveal new physics beyond the SM. Huge data in the 
B-sector has already been accumulated at both the B-factories (Belle and BABAR). This in 
turn has led, for the first time, to the observation of CP violation in the B system, outside 
the Kaon system. In fact, the angle (3 of the unitarity triangle has been measured from the 
time dependent CP asymmetry of the gold plated B® J/tpK s mode by both Belle and 
BABAR, with almost similar values. The current world average of sin 2(3 is [1] 

(sm2(3) b ^ cds = 0.685 ± 0.032 , (1) 

which is consistent with the SM expectation. With the accumulation of more and more data 
the experimental B physics is now all set to enter the unmatched precision era. Unfortunately, 
we have not seen any clear evidence of physics beyond the SM so far, as far as B-physics is 
concerned. 

Already there are some more measurements available at the B-factories, which are not as 
clean as the sin 2(3 measurement in the golden decay mode B — > J/tpK s , but from the pattern 
of deviation observed it appears that these measurements, in the long run with accumulation 
of more data, may reveal the signature of new physics. One of the modes of this kind is the 
decay mode B — > 4>K S) where in the SM one expects to obtain the same value of sin 2(3 from 
its CP asymmetry measurements, as in the case of B — > J/tpK s , with a correction of 0(X 2 ) 
[2]. The basic difference between these two modes is that the golden mode is tree dominated 
{b — > ccs) whereas the decay mode B — > (pK s is penguin dominated (b — > sss). It should 
be reminded here that earlier the deviation between these two measurements was very large 
but with the accumulation of more data the difference has reduced somewhat. The present 
averaged value is [1] 

(sin2/3) 0i , s = 0.47 ±0.19, (2) 

which has about la deviation from the corresponding cc measurements. In future, even 
if the (sm2(3),pKs value stabilizes around the present central value, with error bars reduced 
then it might show the presence of new physics (NP). Moreover, there are other decay modes 
(involving b — > sss transition) where the data show a similar trend. Except the decay mode 
B° — > rj'K , the value of sin 2(3 extracted from all such modes are within la deviation 
from the corresponding cc value [3]. The present average value of B° — > i]'K° mode is 
(sin 2(3) v > K o = 0.48 ± 0.09, which shows about 2.3 a deviation from Eq. (1). 
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Table 1: Experimental values of Branching Ratio (in units of 1CT 6 ), polarization fractions 
and triple product asymmetries for the — > (f)K*° mode. The Belle results for the triple 
product asymmetries are obtained from combined <pK*° and (pK* + data. 

Observables BABAR[4] Belle [5] Average 

B 9.2 ±0.9 ±0.5 lO.Olljjt&l 9.4 ± 0.9 

f L 0.52 ±0.05 ±0.02 0.45 ±0.05 ±0.02 0.49 ± 0.04 

/ ± 0.22 ±0.05 ±0.02 0.30 ±0.06 ±0.02 0.25 ± 0.04 

Al -0.02 ±0.04 ±0.01 §[0.01 ± 0.10 ±0.02] -0.01 ±0.03 

A° T 0.11 ±0.05 ±0.01 J[0.16±S:i4 ± 0.03] 0.10 ±0.04 



The vector- vector counterpart of the seemingly problematic B — > <pKs decay mode, i.e., 
B — > (f)K*, governed by the same b — > sss transition as in B — > (f)Ks, has also created a 
lot of attention recently. Both BABAR [4] and Belle [5] have observed this decay mode and 
the measured quantities are summarized in Table- 1. The measured longitudinal polarization 
fraction in this mode is well below from its expected value [6], i.e., /l ~ 1 — 0(l/ml), 
widely known in the literature as the polarization anomaly in B — > <f)K*. In the B rest frame 
both the vector mesons are emitted back-to-back and from the spin angular momentum 
conservation it follows that both the vector mesons are paired up with the same helicity 

combinations (like 00, ±± and , i.e., the helicity combinations out of the three possible 

helicity states for each vector meson, namely, A = 0, ± and — ). In the SM, it so happens 
that the helicity combination 00, called longitudinal, (i.e., for both the vector mesons the 
spin direction is proportional to the direction of motion ) is almost the only preferred one 
and the occurrence of other two possible helicity combinations is suppressed by 0(l/m 2 B ), 
itlb being the B-meson mass. Thus, the longitudinal polarization fraction is defined as the 
ratio of the decay rate corresponding to the longitudinal polarization (say I\ ) to that of 
the total decay rate T, i.e., fi = T L /T pa 1. However, as seen from Table- 1 its measured 
value is only about ~ 50% of the expected value. 

The unexpected deviation of /& from the expected value of (9(1) is known as the polariza- 
tion anomaly in B — > <pK* decay. In practice, the value of /l is slightly less than unity (and 
predicted to be around 0.9) in the SM. We would like to mention here that the polarization 
measurement in all other vector- vector modes, observed so far, (e.g., B — > pK* and B — > pp) 
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are in accordance with the SM expectations. 

Speculation of the existence of new physics in B — > <pK* was pointed out in [7], in the 
context of two beyond the standard model scenarios (namely, R-parity violating supersym- 
metry and vector like down quark model). The issue of longitudinal polarization problem in 
B — > <f)K* process and its implications were outlined in the review talk [8]. Recently, there 
have been a lot of works in this issue that one can find in the literature [9, 10, 11, 12, 13, 14, 15] 
both in and beyond the standard model. To be more specific, what we essentially need is a 
destructive longitudinal component, or an enhancement in the transverse component or the 
occurrence of both of them coherently to account for the low longitudinal polarization (or 
large transverse polarization) observed in B — > <f)K*. But, at this stage, because of our lack 
of complete understanding of quark-hadron dynamics, it will be very hard to pinpoint the ex- 
act nature of it. Nevertheless, at least, it will be immensely rewarding to see if one can really 
afford to have similar behaviour in some of the popular beyond the SM scenarios. There- 
fore, in this paper we would like to analyze the decay mode B^ — > <pK* in one of the most 
popular beyond the SM scenarios, i.e., in the minimal supersymmetric standard model with 
mass insertion approximation [16, 17] and to see whether the observed polarization anomaly 
can be accounted for in this model or not. It should be noted here that the contributions 
arising from gluonic dipole operator with squark-gluino loop, are enhanced by a factor of 
(rrig/mb) compared to the SM contributions due to the chirality flip from the internal gluino 
propagator in the loop and it will interfere destructively with the longitudinal component 
of SM amplitude for RL mass mass insertion (which is what exactly one needs, as noted 
earlier). Therefore, one would naively expect that the polarization anomaly in B — > <pK* 
mode can possibly be explained by the minimal supersymmetric standard model with RL 
mass insertion. 

The paper is outlined as follows. In the next section we present the basic formalism 
for the B — > V1V2 decay mode. In section 3, we calculate the SM contribution in QCD 
factorization approach for the sake of completeness. Section 4 contains the new physics 
contribution to account for the lower longitudinal polarization and in section 5, we present 
our conclusions. 
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2 Polarization fractions and Triple Product asymme- 
tries in B —> V1V2 decay 

The most general covariant amplitude for the decay mode — > V1V2 can be described as 
[19] 



A(B° d (p) -> ^(pi,ei)V 2 (p2,e 2 ) = 



mim2 mim 2 



(3) 



where p is the 5 meson momentum and rrii, pi and Si (i = 1, 2) denote the masses, momenta 
and polarization vectors of the outgoing vector mesons. 

However, it is customary to express the angular distribution of B® — > V1V2, with each 
vector meson subsequently decaying into two particles, in terms of the helicity amplitudes 
usually defined as 

^A = (K(A)V 2 (A)|7< e// |SS>, (4) 

for A = 0,±1. 

The relationship between the helicity amplitudes and the invariant amplitudes a, b, and 
c are given as 

H±i — a ± c\J x 2 — 1 , H = —ax — b(x 2 — 1) , (5) 

where x — (pi ■ p2)jva\in2 = (m 2 B — mf — ml) / (2miui2)- The corresponding decay rate using 
the helicity basis amplitudes can be given as 

where p c is the magnitude of c.o.m. momentum of the outgoing vector particles. It is also 
convenient to express the relative decay rates with longitudinal and transverse polarizations 

as 
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(7) 

The helcity amplitudes H\ for the decay B® — > V1V2, where V\ and V2 are the antiparticles of 
V\ and V2 respectively, have the same decomposition as (3) with a — > a, b — > b and c — > — c. 
a, 6, and c can be obtained from a, 6 and c by changing the sign of their weak phases. 
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To take advantage of more easily extracting the CP odd and CP even components the 
angular distribution is often written in the transversity basis. The amplitudes in transversity 
and helicity basis are related to each other through the relations 



A = 



H+i — H_i 



A,, = 



H +1 + H_ x 



A = H . 



V2 ' " " V2 
In the transversity basis the longitudinal and the CP-odd polarizations are given as 
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(8) 



(9) 



The triple product asymmetries (TPA's) in — > ViV 2 decays are defined as [20] 

A = 



Ai = 



lm(A ± A* ) lm(A ± A*) 



J2x\A> 



Ea \A> 



]m.(A ± AX) lm(A ± Al 



(10) 



where A = 0, II, ±. 



3 Standard Model contribution 

In the SM, the decay process B° d — >■ (j)K*° receives contribution from the quark level transition 
b — > sss, which is induced by the QCD, electroweak and magnetic penguins. The effective 
Hamiltonian describing the decay b — > sss [21, 22] is given as 



G, r 10 



y/2' 



Heff = —F=V nh V*, 



11) 



£ cp 3 + c g o g 

-j=3 

where q — u, c. O3, • • • , Oq and O7, • • • , O10 are the standard model QCD and electroweak 
penguin operators respectively, and O g is the gluonic magnetic penguin operator. The values 
of the Wilson coefficients at the scale /j, ks m b in the NDR scheme are given in Ref. [23] as 

Ci = 1.082 , C 2 = -0.185 , C 3 = 0.014 , C A = -0.035 

C 5 = 0.009 , C 6 = -0.041 , C 7 = -0.002a , C 8 = 0.054a , 

C 9 = -1.292a , C10 = 0.263a , C g = -0.143 . (12) 

We use the QCD factorization approach [21, 22] to evaluate the hadronic matrix elements, 
which allows us to compute the nonfactorizable corrections in the heavy quark limit. Naive 
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factorization is recovered in the heavy quark limit and to the zeroth order of QCD correc- 
tions. The decay mode B — > <pK* has been analyzed in Refs. [15, 18, 24] using the QCD 
factorization approach. We will first briefly discuss the essential differences between these 
three approaches. It has been shown in [24] that the magnetic dipole penguin will con- 
tribute to all the three helicity amplitudes with almost same order. Later, this result has 
been changed in Ref. [18] where they have shown that the magnetic penguin will contribute 
only to longitudinal polarization amplitude (Hq). However, very recently again it has been 
corrected in Ref. [15] that the positive helicity amplitude (H + i) also receives a small but 
nonzero contributions from the the magnetic dipole operator. It should be noted here that 
the contribution of the magnetic dipole operator to the tranverse amplitudes are also found 
to be small but nonzero in pQCD approach by Li and Mishima [14]. 

Here, we will use the results of QCD factorization method as obtained in Ref. [15]. In 
this approach the helicity amplitudes are given as 



a ., 



m 



K" 



Am 2 B p 2 



rriK* 



H ±1 = 



m B 
G F 

71 v "'^ 



A BK> 



ml)(m B + m K *)A* K *(ml) 



BK* 



, (13) 



Tn B + m K * 
0,±1. Aff{q 2 ) and V BK * (q 2 ) are the 



where a h = a| + a\ + a\ — \{p!] + dg + a^) with h 
form factors describing the B — > K* transitions [25] evaluated at q 2 = m^. The expressions 
for the effective parameters otl appearing in the helicity amplitudes (13) are given as [15] 
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where A g = V^V^*, CV = (iV 2 — 1)/2N and = m\jm\. The QCD penguin loop functions 
G , ' l (s) are given by 

G°(s) = ^-^ln^ + 4/ 1 ^$[(x)^, S ), 
3 3 m b Jo 11 

G^s) = l-hn^ + 2[ 1 dx ( g W {x )±l*£!M)g {x , 8 ), (15) 



3 3 m b Jo \ 4 dx 

with the function g(x, s) defined as 

g(x,s)= / du u(l — u) In [s — u(l — u)(l — x) — ie] . (16) 
Jo 

The EW penguin type diagrams induced by the operators 0\ and O2 are 

c " = (T t Gh(Su) + Y t Gh(Sc) ) (° 2 + §) • (17) 

The gluonic dipole operator O g gives a tree level contribution as 



G ° = - 2 L ix T 



X 



G + g = J dxlg<?+(x) + 



{v)4> ldg^jx) 



4 dx I 1 — x 
G~ = 0. (18) 

The vertex correction factors fj are given as 

f?(a) = -12 In -18 + 6(1- a) + ^ dx ^Ux) U— —\nx-Zin 
m b Jo " \ 1 — x 



f±\a) = -12 In— -18 + 6(1 -a 

(a)4> 



+ Cdx ff f)*( l)± ^ f 3 i^ln,-3^ . (19) 
jo \ 4 ax /VI — x / 

The hard spectator interaction fj T arising from the hard spectator interaction with a hard 
gluon exchange between the vector meson and the spectator quark of the B meson is given 

as 



N ho Jo p Jo v Jo u 



N m B h±i Jo p Jo 



v 2 



X 



I 1 du (g^(u) - a J^(u)\ ^_2if B f K .Um K .m, f *f(p) 
Jo I 4 / JV m|/i±i Jo p 



x / dfdw 



(V) 



dv 



g ^ iu)± a d 9 ^(u) \u + v^ (20) 



d« 



UV " 



with t> = 1 — v and 

if* 



ho = 
h±i = 



2m,K* 



[m B -m K * - mJ(m B + m K *)A 1 (raj 



1 "'/;/':' .^mv.J 



rag + m K * 



m l) 



{m B + m K *)Af K *{ml) =f 



2m B p c 



(21) 



ra B + ra^* 

The asymptotic form of the leading twist ($J(x), &±(x)) and twist-3 (g±\x), g ( f > (x)) light 
cone distribution amplitudes are defined as 



§Y(x) = $v( x ) = gf\x) = 6x(l - x) , 



gZ\x) = ±[l + (2x-l) 2 \. (22) 

The light cone projector for B meson in the heavy quark limit can be expressed as [21] 

ifBm B 



M a = 



(i+ 07fc{*f(O + M£(O}' 



(23) 



where £ is the momentum fraction of the spectator quark in the B meson, v = (1,0,0,0), 
n_ = (1, 0, 0, —1) is the light cone vector. The normalization conditions are given as 



Jo Jo 



For our numerical evaluation we use 



o 4 X B 



(24) 



(25) 



with As = 0.46 GeV, which parametrizes our ignorance of B meson distribution amplitudes. 

It should be noted that the presence of logarithmic and linear infrared divergences in ff/ 
implies that the spectator interaction is dominated by the soft gluon exchanges in the final 
states. To regulate these divergences, a cutoff parameter of order Aq C d/^6, with Aq C d=0.5 
GeV has been used. 
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Table 2: The parameters of the form factors describing B — > K* transitions. 





Mo) 


A 2 (0) 


V(0) 


F(0) 


0.294 


0.246 


0.399 


Cl 


0.656 


1.237 


1.537 


c 2 


0.456 


0.822 


1.123 



For numerical analysis, we use the following input parameters. The quark masses ap- 
pearing in the penguin diagrams are pole masses and we have used the following values (in 
GeV) as m u = = m s = 0, m c = 1.4 and = 4.8. The decay constants used are (in 
GeV) fs = 0.161, fx* = 0.217, = 0.231 and = 0.156. The form factors are evaluated 
in the light-cone sum rule analysis [26] where the q 2 dependence is given as 

F(q 2 ) = F(0) exp[ Cl (g 2 /m|) + c 2 (g 2 /m|) 2 ] , (26) 

with the parameters as given in Table-2. The particle masses and lifetime of meson was 
taken from [27]. For the CKM matrix elements, we have used [27] 

| V cb \ = 0.0413 ±0.0015 , p = 0.20 ±0.09, p = 0.33 ± 0.05 . 

With these input parameters, we obtain the branching ratio in the SM as 

B(B° d -> <f)K* ) = (6.34 ± 0.46) x 10~ 6 , (27) 

and the longitudinal and the CP-odd polarizations as 

h = 0.89 , f ± = 0.05 . (28) 

The triple product asymmetries Ap^ (10) are found to be identically zero. 

4 New Physics Contributions 

We now consider the contribution arising from NP. In general the effective AB = 1, NP 
Hamiltonian relevant for the b — > sss transition is given as 
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/NP 



n eff oc 



'£(cro i + crd i ) + c a o g + c g d s 



(29) 



where Oi (0 g ), are the standard model like QCD (magnetic) penguin operators with current 
structure (sb)y-A(ss)v±A and Cf p , C g p are the new Wilson coefficients. The operators 0\ 
(O g ) are obtained from Oi (O g ) by exchanging L <-> R. As discussed in Ref. [9] the NP 
contributions to the different helicity amplitudes are given as 



A NP (B° d 
A NP (B[ 



-*0\ 



<PK* ,o,| 
* 07T ). 



oc Cf p 



qNP 



(30) 



Thus, in the presence of new physics, the different amplitudes can be given as 



M = 

A i = 



SM 



iTVP 



A 



SM 



4- A v 
0,|| + ^0,|| — ^0, 

A SM _|_ — a sm 



e l0iY (r o i| -f i 



l + e^(r ± + f ± ) 



(31) 



where r\, with (A = 0, ||, 1) are the ratio of NP (arising from (AOj and C g O fl part of the 
Hamiltonian) to SM amplitudes, f\ are the corresponding values arising from the C^O, and 
C g O g part. <fi N is the relative weak phase between the SM and NP amplitudes. For simplicity, 
we have assumed a common weak phase for the C and C contributions and zero strong phase 
between the SM and the NP amplitudes. 

Thus the branching ratio is given as 



B(B» -> <j>K 



*0\ 



B 



SM 



1 + 



2 — I ToTTT^ COS 07V 



(32) 



EaI^! m I 2 EaI^T 

where R\\ j0 = r|| j0 — fp and R± = r_|_ + r± and £> 5M denotes the SM branching ratio. The 
longitudinal and the CP-odd polarizations now read as 

\A$ M \ 2 [1 + R 2 + 2R cos^/v] 



h 
f± 



Ex\Ai M \*[l + Rl + 2R x cos<f> N } 
\A s ± M \ 2 [l + Rl + 2R ± cos<f )N } 



J2x\A s x M mi + Rl + 2R x co S( f )N ] ■ 
Furthermore, in the presence of NP, the Triple Product asymmetries (10) are given as 

2(R± - Rq) sin 



(33) 



a\ = 



Ex\Af M \ 2 [l + Rl + 2R x cos ( p N } ' 

2(R ± -R\\)sm<f> N 
J2 x \Af M mi + Rl + 2R x cos ( p N ] ■ 



(34) 
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We now analyze the decay process B% — > <pK*° in the minimal supersymmetric standard 
model (MSSM) with mass insertion approximation. This decay mode receives supersymmet- 
ric (SUSY) contributions mainly from penguin and box diagrams containing gluino-squark, 
chargino-squark and charged Higgs-top loops. Here, we consider only the gluino contribu- 
tions, because the chargino and charged Higgs loops are expected to be suppressed by the 
small electroweak gauge couplings. However, the gluino mediated FCNC contributions are of 
the order of strong interaction strength, which may exceed the existing limits. Therefore, it 
is customary to rotate the effects, so that the FCNC effects occur in the squark propagators 
rather than in the couplings and to parameterize them in terms of dimensionless parameters. 
Here we work in the usual mass insertion approximation [16, 17], where the flavor mixing 
i — > j in the down-type squarks associated with q B and qA are parametrized by (S^ B )ij, with 
A, B = L, R and i,j as the generation indices. More explicitly (Sf L )ij = (V^ MjV^)ij /m~, 
where Mj is the squared down squark mass matrix and rriq is the average squark mass. Vj 
is the matrix which diagonalizes the down-type quark mass matrix. 

Thus, the new effective AB = 1 Hamiltonian relevant for the B^ — ■> <f)K* process arising 
from new penguin/box diagrams with gluino-squark in the loops is given as 



n 



SUSY 
eff 



V2 



i=3 



o 9 + c» p d 9 



(35) 



where 0{ (O g ) are the QCD (magnetic) penguin operators and the Cf p (C g ) are the new 
Wilson coefficients. The operators Oj are obtained from 0{ by exchanging L <-> R. 

To evaluate the amplitude in the MSSM, we have to first determine the Wilson coefficients 
at the b quark mass scale. At the leading order, in mass insertion approximation, the new 
Wilson coefficients corresponding to each of the operator at the scale /x ~ m ~ Mw are 
given as [17, 28] 



C NP 



cr 
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5 — 



6 — 



AG F V tb V t * s m 

4G F V tb V t * s m 

V2a 2 s 
'AG F V tb V t * s m 2 

V2a 2 s 



AG F V a V,;m\ 



2:-! 



2.3 



-fax) - fax) - fax) - fax) 
-fax) + fax) + fax) + fax) 
Y^ + ^B^-^P^-fax) 
-fax) + fax) + fax) + fax) 
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c. 



NP 



~2G F V tb V t * a m% 



, , , t — , (36) 

'23 Vm 6 7 6 V x VJ 

where x = m~/m|. The loop functions appearing in these expressions can be found in Ref. 
[17]. The corresponding C^ p are obtained from C^ p by interchanging L <-> i?. It should 
be noted that the (5^)23 contribution is enhanced by (m~ g /m b ) compared to that of the SM 
and the LL insertion due to the chirality flip from the internal gluino propagator in the loop. 
Therefore, the magnetic dipole operators in supersymmetric model are found to contribute 
significantly. 

The Wilson coefficients at low energy Cf /p (/i ~ m b ), can be obtained from Cf p (M w ) 
by using the Renormalization Group (RG) equation, as discussed in Ref. [23], as 

C(fi) = U 5 (fi,M w )C(M w ) , (37) 

where C is the 6x1 column vector of the Wilson coefficients and U5(/x, My/) is the five-flavor 
6x6 evolution matrix. In the next-to- leading order (NLO), Us(/i, Mw) is given by 

VM M W ) = (l + £^> j) Uf(„, M W ) (l - ^1 jj , (38) 

where Xjf\fi,M w ) is the leading order (LO) evolution matrix and J denotes the NLO 
corrections to the evolution. The explicit forms of U 5 (/i, M w ) and J are given in Ref. [23]. 

Since the O g contribution to the matrix element is a s order suppressed, we consider only 
leading order RG effects for the coefficient C^ p , which is given as [28] 

Cf p (m b ) ~ -0.15 + 0.70 C^ P {M W ) . (39) 

For the numerical analysis, we fix the SUSY parameter as rrig = m g = 500 GeV, a s (Mw) = 
0.119. The absolute values of the mass insertion parameters (^1^)235 with A,B = (L, R) are 
constrained by the experimental value of B — > X s 7 decay [17]. These constraints are very 
weak for LL and RR mass insertions and the existing limits come only from their definitions 
I (^ll rr) 23 1 < 1- The LR and RL mass insertions are more constrained and for instance 
with rrig ~ m g ~ 500GeV, we have \{^lr,rl)22,\ < 1-6 x 10 -2 . In our analysis, we will use the 
above bounds on the mass insertion parameters, i.e., 

K^l^I < 1 and |(^ L ) 23 | < 1.6 x 10- 2 . (40) 



13 



Now substituting the values of the RG evoluted Wilson coefficients C/ vp (m b )'s in Eq. (14) 
we obtain the corresponding a^'s and hence with Eqs. (8), (13) and (31) the amplitudes. 
Assuming that all the mass insertion parameters (6^)23 have a common weak phase, we 
obtain the new physics parameters arising from the LL(LR) and RR(RL) mass insertions 
as 

(^o)ll = (fo)RR < 0.44 , (r ) L R = (f ) RL < 1.3 , 

(r\\) LL = (f\\) R R < 6.2 x 1(T 2 , (r||) Lfl = (t^rl < 7.0 x 1(T 2 , 

(r ± ) LL = (f ± ) RR < 6.0 x 10~ 2 , (r ± ) LR = (f ± ) RL < 7.2 x 10~ 2 . (41) 

Let us now analyze the variation of CP-odd polarization fraction and branching ratio in 
the presence of new physics. We first consider the contributions arising from the LL and 
RR mass insertions. As seen from Eq. (41), these contributions are quite small and it is 
expected that they cannot accommodate the observed large CP-odd polarization. Now using 
the maximum values of r\ from (41), we plot the CP odd polarization (33) and the 
branching ratio (32) versus the weak phase <pN for three different cases {LL, RR and in the 
presence of both LL and RR contributions). It is seen from Figure-1 that, indeed the large 
CP odd polarization fraction cannot be accommodated in these cases although the observed 
branching ratio can be accommodated with LL or RR mass insertions. Next we consider 
the contributions arising from LR, RL and the simultaneous presence of LR and RL mass 
insertions. As seen from Figure-2, in this case the observed CP-odd polarization fraction 
and the longitudinal polarization fraction (/^) can be accommodated with either LR 
or RL mass insertions. The branching ratio can also be accommodated with these mass 
insertions as seen from Figure-3. 



5 Conclusions 

Observation of unexpectedly smaller longitudinal polarization (and large transverse polar- 
ization) in the penguin dominated B — > <pK* mode poses a serious challenge both to the 
theorists and experimentalists in B-physics. This has in turn ignited the desire of revealing 
the existence of new physics beyond the SM. While at present there is no clear indication of 
any NP but with the accumulation of more data, if the longitudinal polarization stabilizes 
around the present central value, i.e., fi = 0.5, with reduced error bars, then this might be 
the first clear evidence of new physics in the b — > s penguin decay amplitudes. The polar- 
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Figure 1: The CP-odd polarization fraction (f±) and the branching ratio of the process 
B® — > <fiK* with LL and RR mass insertions versus the weak phase <p N (in degree). The 
horizontal solid lines represent the experimental central value and the dashed lines represent 
the la range. 

ization measurements in various vector-vector modes undertaken by the ongoing B-factory 
experiments and the experiments to be performed at BTeV and LHC-b will definitely be 
able to guide us to understand the dynamics and possibly provide us a meaningful answer 
to our quest for the existence of new physics. 

We have employed supersymmetry with mass insertion approximation and shown that 
the low longitudinal polarization in B° d — > <f)K*° can be accommodated in this beyond the 
standard model scenario with either LR or RL mass insertion approximation. 

Since in the B-factory data so far we have only seen some kind of deviation in the b — * sss 
transition measurements it may be worthwhile to continue our effort in this direction and 
check carefully if we can really observe NP in this type of penguin induced transitions. If 
NP is present in the b — > sss transitions and indeed if it is responsible for the observed 
lower longitudinal polarization then we expect to see the same effect of lower longitudinal 
polarization in another charmless vector- vector mode, i.e., B s — > 00, which is governed by 
the same penguin induced b — > sss transition. Already the branching ratio for this mode 
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has been measured by the CDF Collaboration [29, 30] and we are looking forward to the 
polarization measurements in this mode. This in turn, at least, will provide us a clear 
picture of the charmless vector-vector transitions induced by the b — > sss penguins and 
possibly revealing the existence of NP in penguin induced b — > sss transitions. If confirmed, 
the polarization anomaly along with the deviation measured in sin2/5 measurements may 
bring us one step further towards the establishment of NP in the penguin induced b — > sss 
transitions. It is therefore urgently needed to closely examine experimentally all the possible 
charmless vector-vector modes to confirm or rule out the existence of new physics. 
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